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Abstract. Supernova 1006 is the first shell type super- 
nova remnant to show evidence of particle acceleration 
to TeV energies. In the present paper we examine this 
possibility by modeling the observed X-ray non-thermal 
emission in terms of synchrotron radiation from Fermi ac- 
celerated electrons. The predicted synchrotron spectrum 
fits the radio and non-thermal component of the observed 
soft X-ray to hard X-ray emission quite well. These par- 
ticles can produce TeV gamma rays by inverse Compton 
scattering on the microwave radiation and other ambient 
fields, and the derived electron distribution is also used to 
calculate the expected inverse Compton flux. We find that 
if the remnant is characterised by a magnetic field strength 
lower than ~ 7/zG, then the TeV flux can be higher than 
that of the Crab Nebula. About 75% of the TeV emission 
from SN 1006 is expected to be concentrated in the syn- 
chrotron bright NE and SW rims (the "hard aegis") of 
the remnant, which would allow a sensitive search if the 
Atmospheric Imaging Cherenkov Technique is used. 
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1. Introduction 

The type la supernova remnant G327.6+14.6, or SN 1006, 
is the remnant of the explosion which took place in AD 
1006. The distance to this source is between 1.4 and 2.1 
kpc (Green 1988). Its mean expansion rate has been mea- 
sured as 0.44 ± 0.13 arcsec/yr, which implies a present 
expansion speed of 3700 ± 1300 km/s for an average dis- 
tance of 1.8 kpc. The time dependence of the expansion 
(R oc i°- 48 ± 013 ) i s consistent with Sedov expansion, or 
with a forward/reverse shock pair moving into constant- 
density material (Moffett, Goss, & Reynolds 1993). 

The ASCA detection of power law X-ray emission from 
the bright northeastern (NE) and southwestern (SW) rims 



of SN 1006 led Koyama et al. (1995) to infer the existence 
of electrons accelerated by the first order Fermi mecha- 
nism up to energies of ~ 200 TeV. The detection of TeV 
7-rays from this system will prove the existence of such ul- 
trarelativistic electrons (De Jager et al. 1995, Mastichiadis 
1996, Pohl 1996). 

In Section 2 we will compile the non-thermal radio to 
X-ray spectrum, which may be due to synchrotron emis- 
sion from these electrons, and in Section 3 we will com- 
pare the acceleration and loss timescale for electrons ac- 
celerated by the first order Fermi process, and show that 
the resulting maximum electron energy predicts a syn- 
chrotron cutoff frequency which is consistent with obser- 
vations. The best-fit spectrum also leads to the derivation 
of the magnetic field strength £?* associated with the shell. 
In Section 4 we will calculate the inverse Compton (IC) 
7-ray spectrum associated with this best-fit model spec- 
trum from Section 3. The detection of TeV 7-rays from 
SN 1006 will enable the observer to determine the mag- 
netic field strength, and an important parameter associ- 
ated with Fermi acceleration from the observed 7-ray flux. 

2. The observed synchrotron spectrum of SN 1006 

The radio morphology shows two bright arcs towards the 
NE and SW (see e.g. Reynolds & Gilmore 1993), with a 
composite spectrum which is given by 
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and shown in Fig. 1. The individual radio flux measure- 
ments for the total emission from SN 1006 (reproduced 
from Fig. 1 of Reynolds 1996) are also shown for compar- 
ison. 

ROSAT imaging observations of SN 1006 shows that 
the soft X-ray remnant is dominated by thermal [kT = 
0.15 keV) emission towards the SE, whereas the X-ray 
counterparts of the NE and SW radio rims (or "hard 
aegis" ) are characterized by a power law spectrum which 
dominates the total emission from the remnant above ~ 1 
keV (Willingalc et al. 1996). Spatially resolved spectral 
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results by ASCA have shown that these bright rims are 
responsible for ~ 75% of the total X-ray flux from SN 1006 
above 1 keV, and the interior of the remnant also contains 
this non-thermal emission, which dominates the thermal 
component above 2 keV (Koyama et al. 1995). The re- 
sults from earlier non-imaging instruments are therefore 
useful to obtain the non-thermal spectrum above 1 keV 
from the entire remnant. For example: TENMA (Koyama 
et al. 1987), EXOSAT (Jones & Pye 1989), and Ginga 
(Ozaki et al. 1994) obtained well-constrained spectral en- 
ergy indices ranging between —2.0 and —2.3 for energies 
between 1.5 keV and ~ 10 keV. This is confirmed by the 
energy index of -2.0 ± 0.2 of the NE rim found by ASCA 
(Koyama et al. 1995). The EXOSAT ME field of view 
was just large enough to cover the entire remnant, with- 
out significant background contamination, and the energy 
spectrum above 1.5 keV (derived from the given energy 
flux) is given by 

F e = 0.066ej^° keV.cm-V^cV- 1 , (2) 

which is shown together with the radio spectrum in Fig. 
1. 

The change from a —0.56 radio energy index to a <~ — 2 
X-ray energy index above ~ 1 keV is indicative of a spec- 
tral turnover or break around ~ 0.25 keV, where the 
spectral energy index should be around —1 (depending on 
the sharpness of the cutoff). In fact, soft X-ray imaging ob- 
servations of the hard aegis by the EINSTEIN SSS (Becker 
et al. 1980) and EXOSAT LE (Jones & Pye 1989) resulted 
in a consistent spectral energy index of —1.2. Willingale 
et al. (1996) used ROSAT PSPC observations to separate 
the thermal and non-thermal energy fluxes of the hard 
aegis and the rest of the remnant, with the same finding 
as Koyama et al. (1995) that the non-thermal soft X-ray 
flux from the hard aegis contributes to ~ 75% of the total 
non-thermal soft X-ray flux of 2.66 x 10 -10 ergs/cm 2 /s 
(for the 0.1 keV to ~ 2 keV range). Assuming the average 
spectral index of -1.2 found by EINSTEIN and EXOSAT 
(Willingale et al. 1996 obtained a spectral index for the 
SW limb only), but using the ROSAT normalization for 
the total non-thermal flux from SN 1006, the resulting 
energy spectrum for this energy range is therefore 

F e ~ 0.047e lre 1 v 2 keV.cm-V^cV- 1 , (3) 
and is indicated by its spectral index in Fig. 1. 

3. First order Fermi acceleration of electrons in SN 
1006 

The non-thermal synchrotron emission of SN1006 from 
shock accelerated electrons has been discussed by 
Reynolds & Chevalier (1981), Ammosov et al. (1994) and, 
more recently, by Reynolds (1996). Here we apply the 
method given in Mastichiadis (1996 - henceforth M96). 
This allows us to treat the problem in a self-consistent 



manner by calculating the electron distribution function 
at each instant in time from the solution of a time- 
dependent kinetic equation for electrons. The physical 
picture is according to the 'Onion-shell- model' (Bogdan 
& Volk 1983) as applied to time-dependent acceleration 
(Ball & Kirk 1992). Therefore we assume that electrons 
are accelerated in the (parallel) shock wave of the expand- 
ing supernova remnant and at each instance there is a flux 
of relativistic electrons which escape downstream and sub- 
sequently radiate. For these we assume to have a power 
law spectrum in energy, i.e. Q e {E,t) = Q c .o(t)E~ s with 
E < E max (t) and an exponential cut-off for E > E max (t) 
(Webb et al. 1984). For the determination of the normal- 
ization factor we followed M96 (after an original sugges- 
tion by Drury 1992) and set Q e ,o(t) oc R(t) 2 u s h(t) 3 where 
R(t) is the radius of the supernova and u s h(t) is the ve- 
locity of the shock. Furthermore, the maximum energy 
electrons can achieve at a time t depends on whether syn- 
chrotron losses become dominant or not and it is given by 
the minimum of (Lagage & Cessarsky 1983, Webb et al. 
1984, M96 and Reynolds 1996) 

E acc (t) - 5.10- 3 A^f^u^tfBt^ erg (4) 
and 

E loss (t) ~ 2.10-\r 1 / 2 A(r c ) 1 / 2 u sh (t)B~ 1 / 2 erg, (5) 

where w s h is the shock velocity in units of km/sec, B is 
the magnetic field strength in Gauss, t yr is the supernova 
remnant age in years, and A(r c ) = (r c — l)/r c (r c + 1) 
with r c being the compression ratio of the shock connected 
with the spectral index s through the relation s = (r c + 
2)/(r c — 1). Also / is the gyrofactor, i.e. the ratio of the 
particle's mean free path to its gyroradius. The case / = 1 
corresponds to a simple Bohm diffusion. 

The relevant energy losses which were included in the 
kinetic equation were adiabatic due to the expansion of the 
remnant, synchrotron and inverse Compton. The only free 
parameter here is the magnetic field as the photon fields 
on which the electrons lose energy by inverse Compton 
scattering are rather well determined (see section 4 below). 

Adopting standard supernova parameters, i.e. a total 
explosion energy of Wsn = 5.10 50 W5.5o erg, an initial 
shock velocity of u s h,in = 7.10 3 U7 j 3 km/sec and an external 
matter density of p = OAp.4 H-atom/cm 3 (Willingale et 
al. 1996) we solved numerically the electron kinetic equa- 
tion which includes a source term and loss terms as these 
were described above (for more details see M96). The ob- 
tained electron distribution function was folded with the 
synchrotron emissivity to obtain the radiated spectrum as 
this should be observed at a time i now ~ 1, 000 years. This 
way we have made numerical fits to the radio/X-ray data 
by treating as free parameters the strength of the mag- 
netic field downstream, the gyrofactor / and the shock 
compression ratio r c . 
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Fig.l shows the best fit to the radio and X-ray spec- 
tra (shown here as full lines) which was obtained for 
B = 3.5/ 2 / 3 fiG and r c = 3.73 (corresponding to s = 2.1). 
This fit holds for /<30. As / increases above this value 
the spectrum starts to break due to synchrotron losses and 
no satisfactory fit could be found. Therefore X-ray data 
alone might not be able to constrain / significantly. More 
constraints on / can be placed from future TeV observa- 
tions as we will show in the next section. 

We proceed now by giving a simple qualitative pic- 
ture which will help in the understanding of our results: 
According to the standard picture of cosmic ray accel- 
eration in supernova remnants (see, for example, Dorfi 
1991) both Q Gj o and E &cc peak close to the transition of 
the supernova remnant from its free expansion to its Se- 
dov phase which occurs at tsd — 410W / 5 1 g O p^ j^ 3 yr. 
For example, E acc oc t for t < t$d and E acc cx i -1 / 5 for 
t > tsd- Similarly, Q e ,o oc i 2 for t < tsd and Q c ,o t 
for t > tsd- Therefore the radiation observed now should 
be dominated by the particles accelerated at t>t$d- Us- 
ing Eqn. (4) we get that the maximum particle energy at 
t ~ t sd is given by 



i c ,n7Df-l -Itt/1/3 -1/3 1/3 

£acc,sd ^ 1.5 10' Bf L g W 5 50 p OA ' u 7 < 3 erg. 



(6) 



g is a factor (of order unity) that reflects the fact that at 
t = tsd the velocity is not any more u s h,in but somewhat 
smaller since deceleration of the shock should have started 
before that time. This in turn implies that the energy i? acc 
as given by Eqn. (||) is, strictly speaking, an upper limit 
of the maximum energy the electrons can attain during 
acceleration. 

Equation (^|) is the maximum energy electrons can 
attain as long as this is smaller than E\ oss (equation 
||). A comparison of the two expressions for t — 
tgd yields that this is indeed the case as long as 

B< 10- 5 / 1/3 5 1/3 ^5 2/ 5oPo.4 /9 W3f7 9 G. If this relation is sat- 
isfied then synchrotron losses do not become important 
at any epoch, since for t > t$d, E acc oc < -1 / 5 while 
E loss (xt- 3 / 5 . 

^From the above discussion it becomes evident that if 
radiation losses do not become important then the syn- 
chrotron spectrum that the current model predicts can 
be fairly well represented by the synchrotron spectrum 
which is derived from a simple electron distribution func- 
tion (power law times exponential cutoff) of the form 
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Fig. 1. Plot of the best fit synchrotron spectrum (thin dashed 
line) and inverse Compton spectra (thin solid lines indicated 
by the /-values as described in Section 3) calculated from 
the time-dependent electron continuity equation (Section 3). 
The open squares are the radio flux measurements reproduced 
from Reynolds (1996). The thick solid lines represent the ob- 
served X-ray spectra corresponding to Eqns. (^), and (^) as 
discussed in Section 2, and indicated by their spectral indices. 
The thick dashed line is the Crab spectrum given by Djan- 
nati-Atal (1995). 



erg. This in turn gives us an efficiency of cosmic ray elec- 
tron production of 



ao2r 1.04 5 -1.04 Ty -O o 65 p -0.35„0.35 



"7,3 



(8) 



dN B N 



-1.56 



E 



-2.12 



exp 



.Jl(JL)l/2 

Eb £>* 



4rrd 2 ,(7) 



with units in total number of electrons per erg while E is 
in ergs. The best-fit value for the magnetic field is given by 
B* = ^f 2/3 g 2/3 W-^ 9 p- 2 4 /9 u 2 7 ^ pG while E b = £ acc , Sd . 
The area factor 47r<i 2 allows us to calculate photon spectra 
in terms of the flux at Earth if {dN / dE) / Aitd 2 is used. Us- 
ing the best fit value for B* the total energy content in rel- 
ativists electrons is ~ 10 49 f~ 1M g' lm W^ P Q 4 35 u^f 



in rough agreement with the hypothesis that about a few 
percent of the supernova energy available for acceleration 
goes to electrons. 

4. The expected 7-ray spectrum of SN 1006 

Several soft photon fields may contribute to the inverse 
Compton (IC) scattering of relativistic electrons into the 
7-ray region. De Jager (1996) discussed the production of 
7-ray spectral features introduced by the scattering of sev- 
eral thermal photon spectra by the electron cutoff (with 
maximum electron energy 7{,mc 2 ) in SNR W44, which may 
explain the EGRET detection (2EG J1857+0118, Thomp- 
son et al. 1996) of this remnant. The cutoff in the spectrum 
derived in Section 3 will similarly produce 7-ray bumps on 
a eF e plot (as in Fig. 1) at 7-ray energies ~ 3kT^ (in the 
Thomson limit) or ~ j b mc 2 (in the Klein-Nishina limit), 
with Ti = 2.76K representing the CMBR, T 2 - 25K the 
galactic dust contribution, X3 = 3, 000 — 4, 000K repre- 
senting the Population II stars, and T 4 ~ 7, 500K rep- 
resenting the contribution from Population I stars, with 
densities given by Skibo (1993). Arendt (1989) also did not 
detect any significant FIR emission associated with swept 
up dust in SN 1006, so that there would be no contribu- 
tion to soft photons from sources inside or near SN 1006, 
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in contrast with W44 where the emission from swept-up 
dust dominates all soft photon fields. 

We include all soft photon fields, using the full cross 
section for isotropic IC as reviewed by Blumenthal & 
Gould (1970) in the calculation of the 7-ray spectrum. 
The IC spectra calculated from the best-fit electron spec- 
trum (Eqn. 0) for different /-values are also shown in 
Fig. 1, and it is clear that the 7-ray energy flux peaks 
at a few TeV (nearly independent of the value of /), the 
dominant source of this bump being the CMBR photons 
scattered in the Thomson limit. The integral flux above 1 
TeV is given by 

F(> ITeV) = 4.1 x io- 11 /- 1 - 35 cnC 2 s _1 , (9) 

with / constrained to values below ~ 20. The Crab spec- 
trum in the range 0.5 to 10 TeV (Djannati-Atai 1995) is 
also shown for comparison in Fig. 1, and it is clear that 
we may expect SN 1006 to be a stronger source than the 
Crab at 1 TeV if f < 3. It is also clear that the EGRET 
upper limit of SN 1006 shown in Fig. 1 does not constrain 
/ significantly, but TeV observations will provide valuable 
limits. 

Finally we would like to note that the TeV flux as 
given above is about an order of magnitude higher than 
the expected flux from the nuclear component of cosmic 
rays assumed to be accelerated also at the supernova shock 
(Drury et al. 1994). 

5. Summary-Discussion 

The recent X-ray observations of SN 1006 (Koyama et al. 
1995) make a strong case for particle acceleration in su- 
pernova remnants. A necessary consequence of this is that 
SN 1006 can be a source of TeV radiation as the acceler- 
ated electrons will scatter off the ambient soft photons to 
very high energies. In the present paper we calculated the 
expected TeV flux by following the basic principles of first 
order Fermi shock acceleration to fit the radio/X-ray ob- 
servations and consequently we used the derived electron 
distribution function to calculate the TeV emissivity. In 
this case we found that the TeV flux depends sensitively 
on the gyrofactor / of the accelerated electrons and con- 
sequently one can use any future TeV observations of SN 
1006 to put limits on this. 

We find that the EGRET upper limit (derived from 
Thompson et al. 1996) does not constrain / significantly, 
but imaging TeV observations have no difficulties observ- 
ing sources nearly 10 times weaker than the Crab. We may 
therefore be able to detect SN 1006 at TeV energies if / is 
not too large. The NE and SW limbs may be detectable by 
the Atmospheric Imaging Cherenkov Technique (see e.g. 
Weekes et al. 1989) if a search for the NE and SW rims are 
made. However, the flux for each rim would be ~ 38% of 
the flux given in Eqn. (j^) for the total TeV emission from 
SN 1006, but by superimposing the pixels from the TeV 



image corresponding to the NE and SW rims, the total 
flux should be about 75% of the flux given in Eqn. (^) . 
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